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CONFINEMENT  AND  REDISTRIBUTION  OF  CHARGES  AND  CURRENTS 


ON  A  SURFACE  BY  EXTERNAL  FIELDS 


ABSTRACT 


The  old  problem  of  light  scattering  from  a  perfectly  conducting  surface 
is  addressed.  An  electromagnetic  field  is  incident  upon  the  boundary,  where 
it  induces  a  charge  and  current  distribution.  These  charges  and  currents 
emit  the  reflected  fields.  A  set  of  equations  for  the  charges  and  currents 
on  the  surface  is  derived  by  eliminating  the  E-  and  B-fields  from  Maxwell's 
equations  with  the  aid  of  the  appropriate  boundary  conditions.  An  explicit 
and  general  solution  is  achieved,  which  reveals  the  confinement  and 
redistribution  of  the  charge  and  the  current  on  the  surface  by  the  external 
field.  Expressions  are  obtained  for  the  surface  resolvents,  or  the 
redistribution  matrices,  which  represent  the  surface  geometry.  Action  of  a 
surface  resolvent  on  the  incident  field  evaluated  at  the  surface,  then 
yields  the  charge  and  current  distributions.  The  Faraday  induction  appears 
as  an  additional  contribution  to  the  charge  density.  Subsequently,  the 
reflected  fields  are  expanded  in  spherical  waves,  which  have  the  surface- 
multipole  moments  as  a  source.  Explicit  expressions  are  presented  for  the 
surface-multipole  moments,  and  it  is  pointed  out  that  charge  conservation  on 
the  surface  sets  constraints  on  these  moments.  The  results  apply  to 
arbitrarily-shaped  surfaces  and  to  any  incident  field.  For  a  specific 
choice  of  the  surface  structure  and  the  external  field,  the  solutions  for 
the  charge,  the  current  and  the  reflected  fields  are  amenable  to  numerical 


evaluation. 


I. 


INTRODUCTION 


The  study  of  chemistry  and  physics  near  a  surface  has  developed  rapidly 

during  the  last  decade.  Investigations  range  from  classical  processes  like 

X  2-4  5  22 

periodic  deposition,  image  formation*  and  dispersion  of  plasmon  waves'3 

to  quantum  mechanical  issues  as  Raman  scattering  of  intense  laser 

23  24  23  26 

light,  *  atomic  fluorescence  near  a  rough  surface,  ’  the  coupling  of 

27 

an  atomic  dipole  to  surface  polaritons  and  cooperative  emission  processes 
28 

near  a  conductor.  It  appears,  however,  that  besides  these  well- 

established  theories,  even  the  simplest  problem  --light  scattering  from  an 

arbitrarily  shaped  surface--  is  not  yet  completely  tractable.  Early 

approximations  like  the  Rayleigh-Fano  expansion  (neglection  of  multiple 

reflections)  or  the  small-roughness  limit  provide  sufficient  understanding 

of  the  induced  effects  on  a  boundary  by  incident  fields,  but  exact  solutions 

in  the  form  of  general  expressions  for  the  scattered  fields  and  the  surface 

waves  are  not  available  at  present.  Contemporary  closed-form  solutions 

pertain  only  to  polarized  plane  waves,  incident  upon  gratings  with  well- 

defined  geometries,  like  square  or  sinusoidal  wells.  The  results  always 

rely  on  the  periodicity  of  the  surface  roughness,  which  implies  the 

applicability  of  Fourier-series  expansions,  or  a  numerical  solution  of  the 

7  8  11 

extinction  theorems,  as  they  exist  in  many  phrasings.  *  *  In  this  paper 
we  consider  a  metallic  surface,  which  is  illuminated  by  an  externally- 
applied  electromagnetic  field  with  an  arbitrary  time-dependence  and  spatial 
distribution.  The  surface  is  not  assumed  to  be  periodic,  and  our  results 
apply  equally  well  to  a  closed  surface  or  to  assemblies  of  surfaces,  as  for 
example  a  sphere  near  a  grating.  We  achieve  closed-form  solutions  of 
Maxwell's  equations  for  the  charge  and  current  distributions  on  the  surface 
and  for  the  reflected  fields,  although  at  the  expense  of  the  assumption  that 
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the  metal  has  a  perfect  conductivity.  The  structure  of  the  expressions  is, 
however,  quite  transparent,  and  it  might  be  expected  that  a  proper  account 
of  the  finite  conductivity  yields  similar  results. 

II.  THE  FIELD  EQUATIONS 

The  time-development  and  the  spatial  distribution  of  the  charge  density 
p(r,t),  the  current  density  j(r,t),  the  electric  field  E(r,t)  and  the 
magnetic  field  B(r,t)  are  governed  by  Maxwell's  equations.  If  we  adopt  a 


Fourier  transform  of  the  real  fields 


lr,t)  *  -Rejdci)fi(r,w)e  iwt. 


(2.1) 


and  similarly  for  the  other  three  fields,  then  the  field  equations  read 
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V-[e(r)E(r)]  ■  p(r), 


V-B(r)  *  0, 


VxE(r)  -  iwB(r)  *  0, 


Vx[p(r)  B( r ) ]  +  iue(r)E(r)  *  j(r), 


(2.2) 


(2.3) 


(2.4) 


(2.5) 


where  we  have  simplified  the  notation  by  writing  E(r)  rather  than  fc(r,io). 
The  frequency  dependence  of  the  fields  and  of  e(r)  and  p(r)  will  be 
suppressed  throughout  this  paper. 

We  shall  suppose  that  the  entire  space  is  occupied  by  two  kinds  of 
media,  perfect  conductors  and  perfect  insulators,  which  are  separated  by 
boundaries.  The  set  of  all  boundaries  will  then  be  referred  to  as  the 
surface.  Within  each  medium  the  dielectric  constant  e(r)  and  the 
permeability  y(r)  will  be  assumed  to  be  r- independent ,  but  across  the 


surface  c(r)  and  p(r)  are  discontinuous.  Conductors  are  specified  by  a 


relation  like  j(r)  *  yE(r),  y>0,  and  the  assumption  of  perfect  conductivity 
implies  the  limit  y-*».  Since  the  current  density  j(r)  should  remain  finite, 
we  obtain  E(r)  =  0  everywhere  in  the  conductor.  From  Eq.  (2.2)  we  then  find 
p(r)  *  0,  and  Eq.  (2.4)  yields  B(r)  =  0,  under  the  restriction  cd  *  0.  In 
this  paper  we  will  exclude  the  trivial  static  case  id  *  0.  Finally,  Eq. 

(2.5)  gives  j(r)  *  0,  and  hence  Maxwell's  equations  in  the  conductor  reduce 
to 

E(r)  -  0,  B(r)  -  0,  p(r)  *  0.  j(r)  -  0.  (2.6) 

Around  a  point  r  on  the  surface  the  fields  are  discontinuous. 
Application  of  Gauss'  theorem  on  (2.2)  and  (2.3)  and  of  Stokes'  theorem  on 
(2.4)  and  (2.5)  enables  us  to  rewrite  the  equations  in  the  vicinity  of  the 
surface  as 

E(r+)  *  e_1o(r)n(r),  (2.7) 

B(r+)  =  pi(r)xn(r).  (2.8) 

Here  o(r)  and  i(r)  are  the  surface  charge  and  current  density,  respectively, 
and  n(r)  represents  the  unit  normal  vector  in  r  on  the  surface,  with  the 
convention  that  it  points  from  the  conductor  to  the  dielectric.  We  have 
introduced  the  notation  r+  to  indicate  a  point  in  the  dielectric  and  close 
to  r.  Explicitly,  we  write 

r+  =  r  +  n(r)6  with  6+0.  (2.9) 

We  note  that  Eqs.  (2.7)  and  (2.8)  combine  the  four  Maxwell  equations  in  r  on 
the  surface,  and  that  they  contain  four  unknown  fields. 

The  dielectric  is  presumed  to  exhibit  no  conductivity  at  all,  so  it  can 
be  specified  by  j  *  yE  with  y-O.  This  implies  j  =  0,  and  from  charge 
conservation  (V‘j  =  iwp)  we  find  p  =  0,  since  we  required  <d  *  0.  Hence,  all 
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charges  and  currents,  if  any,  are  situated  on  the  surface  as  o(r)  and  i(r). 
The  electric  and  magnetic  fields  in  the  dielectric  are  generated  by  o(r)  and 
i(r),  and  they  contain  the  incident  fields.  This  notion  allows  us  to  write 


Maxwell's  equations  for  a  point  r  in  the  dielectric  as 
p(r)  *  0, 


(2.10) 


j(r)  -  0, 


(2.11) 


E(r)  -  E(r)inc  -  ^-JdA'o(r '  )7G(r,r '  )  +  ^JdA'i(r'  )G(r,r' ),  (2.12) 

B(r)  -  B(r)inc  =  ^JdA'i(r* )xVG(r,r' ),  (2.13) 


where  the  integrals  run  over  the  complete  surface.  This  representation 
involves  the  Green's  function  of  the  wave  equation 


G(r.r')  =  | r-r ' | ”  exp( ik|r-r ' | ) , 


(2.14) 


and  its  gradient 

VG(r,r ' )  *  (r-r')|r-r'|"3(ik|r-r' j-l)exp(ik|r-r'|),  (2.15) 

l/2 

which  contain  the  wave  number  k  =  (ey)  w.  We  have  to  solve  the  set  (2.12) 
and  (2.13)  for  o(r),  i(r),  E(r)  and  B(r),  and  Maxwell's  equations  (2.7)  and 
(2.8)  on  the  surface  can  be  considered  as  the  boundary  conditions. 

III.  ELIMINATION  OF  THE  FIELDS 

Maxwell's  equations  in  the  dielectric  medium  are  basically  two 

equations  with  four  unknown  fields,  but  we  can  eliminate  the  radiation 

fields  E(r)  and  B(r)  with  the  boundary  conditions  (2.7)  and  (2.8).  To  this 

end  we  take  r  in  (2.12)  and  (2.13)  as  r+  from  (2.9),  and  then  substitute  the 

boundary  values  for  E(r+)  and  B(r+).  This  procedure  leaves  us  with  a  set  of 

two  equations  for  o(r)  and  i(r).  The  appearance  of  G(r+,r')  and  7  +G(r  »r  ’  ) 

r 
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in  the  integrands  of  (2.12)  and  (2.13)  is  not  convenient  since  it  involves 

points  r+,  which  are  not  situated  on  the  surface.  It  will  turn  out  to  be 

•more  practical  to  have  equations  in  which  the  Green's  function  connects  only 

points  of  the  surface,  rather  than  a  point  on  the  surface  to  a  point  in  the 

dielectric.  However,  care  should  be  exercised  in  replacing  r+  by  r,  because 

the  integrals  are  discontinuous  across  the  surface.  If  we  take  the  limit 
r+->x  properly  (see  Fig.  1),  we  obtain,  for  instance. 


jdA'  i(r '  )xVG(r+,r '  )  *  2iri(r)xn(r)  +  jdA'  i(r '  )x7G(r,r '  ), 


(3.1) 


and  we  observe  that  replacing  r  by  r  requires  that  we  should  add  the  term 

29 

2ni(r)xn(r) .  It  was  already  pointed  out  by  Maradudin  that  integrals  of 

this  kind  appear  to  have  a  finite  contribution  from  a  single  point.  This 

feature  can  however  also  be  regarded  as  resulting  from  the  discontinuity  of 

the  fields  across  the  surface.  Critical  comments  on  this  issue  have  also 

12 

been  made  by  Agarwal  in  a  slightly  different  context. 

With  the  explicit  forms  of  G(r,r')  and  its  gradient,  as  they  are  given 
by  (2.14)  and  (2.15),  we  can  evaluate  the  limit  r+-*-r,  which  yields  the  set 
of  equations 

o(r)n(r)  *  ■^dk’oix '  )7G(r,  r '  )  +  i  (r  '  )G(r,r  ' ) 


+  2E(r)inC, 


(3.2) 


i(r)xn(r)  =  ^“JdA' i(r' )xVG(r,r ' )  +  2y  *B(r)*nC, 


(3.3) 


for  o(r)  and  i(r).  We  can  write  o(r)  and  i(r)  in  the  integrands  as 


o(r)  =  n( r ) • (o( r )n( r ) ) , 


(3.4) 


i(r)  =  n(r)x(i(r)xn(r)). 


(3.5) 
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since  i(r)  is  parallel  to  the  surface,  which  shows  that  Eqs.  (3.2)  and  (3.3) 
are  essentially  a  set  of  equations  for  the  vector  fields  o(r)n(r)  and 
i(r)xn(r)  on  the  surface. 

Equation  (3.3)  for  i(r)xn(r)  has  the  form  of  an  inhomogeneous  Fredholm 
equation  of  the  second  kind,  where  the  external  field  2y  *B(r)inc  is  the 
inhomogenety.  In  the  same  fashion,  Eq.  (3.2)  has  2eE(r)inc  (and  the  current 
term)  as  an  inhomogeneous  part.  Hence  the  incident  fields  can  be  regarded 
as  the  source  terms  of  these  equations.  In  this  sense  a( r)  *  0  and  i(r)  * 
are  a  result  of  the  presence  of  the  driving  field,  so  the  charges  and  the 

currents  are  confined  on  the  surface  by  the  field.  If  there  is  a  net  charge 

on  the  surface,  this  mechanism  might  also  be  conceived  as  a  redistribution 
process.  Equations  (3.2)  and  (3.3)  resemble  the  extinction  theorem  for  the 
analogous  problem  of  scattering  of  an  incident  field  from  a  dielectric 
grating.  The  extinction  theorem  is,  however,  a  homogeneous  equation,  and 
its  solvability  condition  is  equivalent  to  the  dispersion  relation  for 
surface  polaritons. 

IV.  REPRESENTATION  OF  THE  SURFACE 

Ordinary  Fredholm  equations  are  single-variable  equations  for  a 
function  on  the  complex  plane,  and  they  can  be  solved  by  an  expansion  of  the 

function  onto  a  suitable  complete  set.  Our  equations  for  o(r)n(r)  and 

i(r)xn(r)  are  three-dimensional  and  surface-related  equations  for  a  vector 
field,  so  we  have  to  modify  the  standard  techniques  slightly.  In  order  to 
accomplish  this,  we  introduce  spherical  coordinates  (r,6,$)  with  respect  to 
an  arbitrary  origin,  and  we  will  abbreviate  the  direction  6 , d>  by  the  single 
variable  Q.  Then  the  assembly  of  all  points  r,  which  constitute  the 
surface,  can  be  represented  by  a  set  of  functions  £(n)^.  The  £(0) ^  will 
indicate  the  distance  from  the  origin  to  a  point  r  on  the  surface,  in  the 


o  i 


direction  Q,  while  the  subscript  A  accounts  for  the  multiplicity  (see  Fig. 
2).  In  this  fashion,  the  surface  is  divided  in  regions,  numbered  by  A, 
where  its  shape  is  defined  by  a  function  which  determines  uniquely 

the  spherical  coordinates  (r,8,$)  ■  (£(0  , $) A , 8 ,$)  of  a  point  r  in  this 
region.  The  shape-functions  £(£l)A  will  be  assumed  to  be  given,  and 
therefore  we  can  represent  a  point  on  the  surface  by  its  surface-coordinates 
(A,n)  rather  than  by  its  spherical  coordinates  (r,n).  We  will  use  A  as  a 
subscript  and  ft  as  a  variable. 

The  measure  dA(fl)^  and  the  direction  n(fl)A  of  the  surface  at  a  given 
point  (A,fl)  are  fixed  by  its  shape  £(fi)A<  For  instance,  the  infinitesimal 
surface  area  at  ( A , £7 )  is  given  by 

dA(£l)A  =  f(n)xdfl  (4.1) 

with 

f(n)A  -  e(n)AU(n)A  +  (ff5(n)x)2  +  — S- (afc(n)A)2>1/2*  (4,2) 

A  sin^0  v 

in  terras  of  the  infinitesimal  surface  area  d£l  *  sin0d0d$  of  the  unit  sphere. 
Hence  the  function  f(fl)  accounts  for  the  deviation  of  the  surface  curvature 

A 

from  the  curvature  of  a  sphere,  and  with  the  aid  of  (4.1)  we  can  transform  a 
surf ace- integral  over  the  region  A  into  an  integration  over  a  part  of  the 
unit  sphere.  We  note  that  not  every  direction  fl  for  a  given  A  corresponds 
to  a  point  on  the  surface.  It  will  turn  out  to  be  convenient  to  extend  the 
definition  (4.1)  of  f(fi)A  as 

f(0)A  =  0  if  [J  does  not  correspond  to  a  point  on 

the  surface  in  region  A.  (4.2) 

Then  we  can  write  the  surface  integrals  as 


9 


JdA. .  . 


u 


dnf(n)x 


(4.3) 


where  the  integrals  now  run  over  the  complete  unit  sphere  for  every  A.  This 
construction  will  enable  us  to  apply  the  general  theory  of  expanding  vector 
fields  on  a  sphere. 


V.  EXPANSION  OF  THE  FIELDS 

Since  we  are  using  spherical  coordinates,  the  spherical  harmonics 

Y(0)^  supply  a  suitable  complete  set  on  the  unit  sphere  for  an  expansion  of 

the  magnitude  of  a  vector  field.  The  direction  of  a  vector  will  be  expanded 

onto  a  space-fixed  set  of  three  unit  vectors,  denoted  by  e^,  which  is,  for 

instance,  the  Cartesian  set  e  ,  e  ,  e  or  the  spherical  set  e.,,  e„,  e 

-x  -y  -z  ~+l  -0  --1 

Then  the  vector  fields  e  constitute  a  complete  set  on  the  unit  sphere 

tm-x 

for  an  expansion  of  an  arbitrary  vector  field. 

It  is  our  aim  to  solve  Eqs.  (3.2)  and  (3.3)  for  o(r)n(r)  and  i(r)xn(r). 
We  thus  start  with  an  expansion  of  these  fields. 


f(n)Ao(n)An(n)A 


■[ 

£mx 


S„  Y  (n)e  , 
£mxA  £m  -x 


(5.1) 


f(n).  i(n).xn(n).  =  /,  I.  .Y.  (n)e  ,  (5.2) 

A-  A  -  A  £mxA  im  -x 

£mx 

and  note  that  we  have  included  a  factor  f(^)^  on  the  left-hand  side.  This 

is  necessary,  since  otherwise  the  left-hand  side  of  Eqs.  (5.1)  and  (5.2) 

would  not  be  properly  defined  for  every  fi.  The  driving,  incident  fields 
inf'  ^ 

E(r)  ,  B(r)  in  Eqs.  (3.2)  and  (3.3)  enter  only  through  their  value  on 

the  surface,  so  that  we  can  expand  them  on  the  surface  set  according  to 


f(n)AE(n)Ainc 

f(n)AB(n)Anc 


A 

£m 

tmx 


E .  ,Y(n)„  e  , 

£mxA  Urn-x 


B  .Y(n)0  e  . 
UrnxA  £m-x 


(5.3) 


(5.4) 


0 


The  expansion  coefficients  for  the  incident  fields  then  follow  from  the 
inverse  relation 


-  Jdnf(fi)A!(n)AnC-?xY(n)L  (5-5) 

"tfljxA  ■  (5.6) 


and  the  appearance  of  f(ft)^  in  the  integrands  reflects  that  we  actually  have 
integrals  over  region  A  of  the  surface.  This  illustrates  that  f (n)^Y(fl).  e* 
can  be  considered  as  a  complete  surface  set  for  the  expansion  of  a  vector 
field  on  the  surface.  Note  that  we  allow  e^  to  be  complex,  which  is  the 
case  for  a  spherical  set. 


VI.  THE  CHANGE  AND  THE  CURRENT  DISTRIBUTIONS 

It  is  straightforward  to  rewrite  Eqs.  (3.2)  and  (3.3)  for  o(r)n(r)  and 
i(r)xn(r)  in  terms  of  their  expansion  coefficients.  We  obtain 


l 


(R 


,  -6  „ .  ,  6  M51W)S„,  . 


. 'm' t ' A ' 
iujeu 


(2) 

tmxA ,  l 'm' t  '  A '  “ ll '  '"mm  ’^AA  *  "'"l 'm'  x  '  A  '  * 
r 

(3) 


l 


l 'm' x ' A 

,(D 


R  I 

,  tmxA.fc'm'x'A'  H'm'x'A' 


2eE 


Imx  A  ’ 


(6.1) 


t 'm' x ' A 


(R  V  ,  „ .  .  , , ,-600 ,6  ,6  ,6)1,  .  . 

tmxA, l  m  x  A  tt  ran  xx  AA  l  m  x  A 


tmxA  ’ 


(6.2) 


which  are  two  coupled  inhomogeneous  linear  equations  for  the  surface  charge 


and  the  surface  current  I„  The  expansion  coefficients  E„ 

0  mr  \  *  0 


and 


tmxA - -  ‘tmxA'  - . —  “tmxA 

^tmxA  ^or  t*ie  externa^  fields  are  supposed  to  be  given.  The  set  (6.1)  and 

(6.2)  also  involves  three  R-matrices,  with  matrix  elements 


.  ...  =  ^•JdnJdn'f(n).Y(n)*  Y(n 

tmxA,  l  m  x  A  2vJ  J  A  tm 


') 


t '  m ' 


{e^x(n(n'  )A-xe^  ,  )}  -VG(fi,n'  )^.  , 


(6.3) 


v-iv; 
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r(2) 

t' A' 


tiWdQ' 


(n(n')A<.e*Je*.VG(n,fi')AA., 


(6.4) 


C,«v,T  ■  sJdnJ^'f(nV(n)^v(n-)t.n. 


eT-(n(n')A,xeT,)G(n,fl')AA., 


(6.5) 


where  we  have  written  G(n,fl')AA«  for  the  Green's  function,  which  connects 
the  points  (A,fl)  and  (A', O')  of  the  surface.  We  emphasize  that  these  R- 
matrices  depend  only  on  the  geometry  of  the  surface,  and  not  on  the  external 
fields.  Prescription  of  the  shape  of  the  surface  determines  the  R-matrices. 
Recall,  however,  that  the  R-matrices  depend  on  the  frequency  u>  through  the 
Green's  function,  but  this  is  merely  a  parametric  dependence  and  independent 
of  the  external  field. 

The  expansion  coefficients  SjjnxA  can  always  be  arranged  in  a  one¬ 
dimensional  array,  considered  as  a  vector,  and  similarly  R^\  R^2^  and  R^3^ 
can  be  regarded  as  two-dimensional  matrices.  Then  we  can  write  (6.1)  and 
(6.2)  as 


(R(2)-1)S  -  i(D€pR(3)I  -  2eE, 
(R(1)-1)I  -  -2p'1B, 


(6.6) 

(6.7) 


where  we  have  also  adopted  a  vector  representation  for  the  driving  fields. 
The  solution  of  (6.6)  and  (6.7)  is  immediately  found  to  be 


s  =  ■  iujR(3) — nrB}’ 

1-RV  '  1-RV  ' 


•1 


1-R 


(1) 


B, 


(6.8) 

(6.9) 
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which  expresses  the  charge  density  S  and  the  current  density  I  explicitly  in 
the  externally  applied  fields  E  and  B  and  the  surface-shape  matrices 
R(2)  and  R(3\ 

For  vanishing  external  fields,  e.g.,  E  *  0  and  B  ■  0,  the  charge  and 
current  distributions  also  vanish,  as  can  be  seen  explicitly  from  Eqs.  (6.8) 
and  (6.9).  Hence  the  charges  and  currents  are  indeed  confined  to  the 
surface  by  the  external  fields.  Remember  that  we  have  excluded  the  static 


case  <d  ■  0,  for  which  we  can  have  charges  on  a  surface  without  external 

(2)  -1 

fields.  Furthermore,  we  can  identify  the  resolvents  (1-R  )  and  (1- 

R^^)  3  as  the  operators  which  account  for  the  redistribution  of  the  charges 

and  currents,  respectively,  as  resulting  from  the  Lorentz  force  between 

charges  and  between  currents.  The  coupling  of  charges  and  currents,  which 

(3) 

is  the  Faraday  induction,  is  incorporated  in  the  R  -matrix. 


VII.  THE  REFLECTED  FIELDS 

The  incident  field  induces  charges  and  currents  on  the  surface,  and 

these  oscillating  charges  and  currents  emit  radiation,  which  are  the 

reflected  fields.  In  this  section  we  express  these  fields  in  terms  of  the 

expansion  coefficients  S„  ,  and  I.  . ,  as  they  are  given  explicitly  in  the 

unxA  iraTA 

previous  section. 

In  Eq.  (2.12)  we  expressed  the  reflected  electric  field  E(r)-E(r)inC  in 
terms  of  o(r)  and  i(r),  and  similarly  B(r)-B(r)*nC  in  terms  of  i(r).  With 
(3.4)  and  (3.5)  we  can  rewrite  these  equations  in  a  way  that  c(r)n(r)  and 
i(r)xn(r)  are  the  source  fields,  and  then  we  can  apply  (5.1)  and  (5.2)  in 
order  to  find  an  expansion  on  the  spherical  set.  However,  the  resulting 
expressions  are  not  transparent,  since  they  will  involve  the  Green's 
function  and  its  gradient.  In  order  to  achieve  a  more  comprehensible 

23 

result,  we  expand  the  Green's  function  on  the  spherical  set.  We  write 
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[r.r')  -  4irik^  h(l)(k£(fi')  )  Y(n')  j(kr)  Y(fl)* 


(7.1) 


where  and  j  are  spherical  Bessel  functions.  Here  the  convention  is 

that  we  choose  the  origin  of  our  coordinate  system  in  the  dielectric,  and  in 
such  a  way  that  the  inequality 


c(n)A  >  r 


(7.2) 


holds  for  every  (A,0).  The  vector  r  is  the  position  in  the  dielectric, 
where  we  wish  to  evaluate  the  reflected  fields.  The  expansion  coefficients 
^tmxA  an<*  ^£mxA  ^ePen£*  on  t*ie  position  of  the  origin,  so  both  the  charge  and 
current  distributions  and  the  reflected  fields  must  be  evaluated  with 
respect  to  the  same  coordinate  system.  Furthermore,  restriction  (7.2)  must 
hold  in  order  to  apply  the  series  expansion  (7.1)  of  the  Green's  function. 
For  a  given  r,  this  can  always  be  arranged. 

The  solution  for  the  fields  can  be  cast  in  an  appealing  form  by  the 
introduction  of  the  source-term  vectors 


s|»  -  Is,  .e 

~£m  u  £mxA-x 
x 

I(A)  -  ll  e  . 
=lm  L  £mxA~x 


(7.3) 


(7.4) 


In  view  of  (5.1)  and  (5.2),  these  and  1^^  are  just  the  expansion 

coefficients  of  f (O)^o(n)^n(n)^  and  f (n)Ai(fl)Axn(n)^  after  an  expansion  of 

these  fields  onto  the  set  of  spherical  harmonics,  but  without  a 
decomposition  along  the  basis  vectors  e^ .  Furthermore,  we  define  the  vector 


-  £m ,  £ '  m  * 


|dnf(n)AY(n)imY(n)i>m,hv  '(k5(n))rn(n)A, 


(7.5) 


which  is  a  surface  integral  over  the  region  A.  It  is  the  integrated  normal 
vector  n(fi)^  times  the  appropriate  weight-functions.  This  vector 
depends  only  on  the  shape  of  the  surface.  After  these  preliminary 


definitions,  we  can  write  for  the  reflected  fields 


!r)-E(r)ino  -  4 


x<m  x.  m  A 


"iU)  j  ».  j(krVY(nVm'  ’ 

Mu  A  ’ 


(7.6) 


B(r)-B(r): 


tat  m  A 


(7.7) 


These  explicit  expressions  for  the  fields  that  are  emitted  by  the  surface 

charge  and  current  distributions  exhibit  a  clear  separation  between  the 

source  terms  s!|^  and  if^and  the  redistribution,  due  to  the  surface 
-£m  ~fcm 

geometry,  which  is  accounted  for  by  the  vector  P»*\.  ..  The  spatial 

distribution  is  represented  as  an  expansion  in  the  spherical  waves 
jOcr^YOJ)*^  and  7j (kr) ^(fl)*. 

VIII.  SURFACE  MULT I POLES 

We  can  eluciate  the  significance  of  the  expansions  (7.6)  and  (7.7)  for 
the  reflected  fields  by  the  introduction  of  the  surface  multipoles.  To  this 
end  we  define  the  multipolar  moments  of  the  charge  and  the  current 
distributions  as 


C  -  W  D(X)  .<5(X) 

lm  e“,  m ', £m  -  £ '  m '  ’ 

i  m  A 


*“  tfc., 


)  PM  xl(*> 

Zl'm'  ,lm  -I’m' * 

i  ffl  A 


(8.1) 


(8.2) 


where  C„  is  a  scalar  and  J„  is  a  vector.  These  multipolar  moments 
£ra  ~£m 

represent  the  charge  and  current  distribution  of  the  complete  surface,  not 


just  in  one  region  A.  The  emitted  fields  now  attain  the  form 


E(r)-E(r)inc  -  \ 


L  Ctm7J(kr)£Y(n)tm  '  -£m 

£m  £ra 


l 


V(n)ta- 


(8.3) 
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B(r)-B(r)inc  -  [  J  x7j(kr)  YOl)^.  (8.4) 

which  greatly  resembles  the  multipole  expansion  of  the  fields  emitted  by  a 

charge  and  current  distribution  in  a  restricted  region  of  space.  The 

distinction  is  of  course  that  the  source  terms  C»  and  J.  here  gain 

im  -im  * 

contributions  from  everywhere  in  space,  rather  than  from  a  localized  area. 

This  results  effectively  in  an  exchange  of  the  spherical  Bessel  function 

h^^(kr),  with  j(kr)^  in  the  expansion  of  the  Green's  function. 

The  surface  multipolar  moments  C.  and  J.  are  not  independent.  Tron* 

un  -Jun 

the  fact  that  the  fields  obey  Maxwell's  equations,  as  they  do  by 

construction,  it  follows  that  they  are  subject  to  some  constraints. 

1.  DC 

From  V*(E(r)-E(r)  )  *  0  we  readily  derive  the  relation 


+  £±E[ 

V2I+3  p=*(i+l)  T 


(8.5) 


for  a  spherical  basis  set  e^.  Here  ( imlx | itlp)  denotes  a  Clebsch-Gordan 
coefficient.  The  constraint  (8.5)  can  be  considered  as  the  surface- 
integrated  form  of  charge  conservation  (V  -j  =  iwp)  for  the  surface  charge 
density  o(r). 


IX.  CONCLUSIONS 

We  have  studied  the  charge  and  current  distributions  on  the  boundary  of 
a  perfect  conductor  with  a  dielectric,  as  they  are  cconfined  and 
redistributed  there  by  an  externally  applied  electromagnetic  field.  The 
surface  was  allowed  to  have  an  arbitrary  shape,  and  we  did  not  impose  any 


periodicity  condition.  We  obtained  closed-form  and  exact  expressions  for 
o(r)  and  i(r)  everywhere  on  the  surface.  This  was  accomplished  by  deriving 
a  set  of  inhomogeneous  Fredholm  equations  of  the  second  kind  for  o(r)n(r) 
and  i(r)xn(r)  from  Maxwell's  equations,  and  subsequently  solving  these 
equations  by  an  expansion  on  a  discrete  spherical  set  of  basis  vector 
functions.  The  solution  involves  surface-structure  matrices,  the  R- 
matrices,  which  are  independent  of  the  incident  field.  It  appears  that  an 
operation  of  a  resolvent  (1-R)  *  on  the  vector  representation  of  the 
impinging  field  on  the  surface  yields  the  charge  and  current  distributions. 
The  Faraday  induction  between  the  E-  and  the  B-fields  gives  rise  to  a 
coupling  between  the  equations  for  o(r)n(r)  and  i(r)xn(r),  and  it  was 
accounted  for  by  the  matrix 

Next,  the  structure  of  the  fields,  which  are  emitted  by  the  oscillating 

charges  and  currents,  was  examined.  The  solution  was  cast  in  the  form  of  a 

spherical  multipolar  expansion,  and  the  multipolar  moments  were  identified 

explicitly  in  terms  of  the  solutions  for  o(r)  and  i(r).  The  effect  of  the 

surface  geometry  could  be  incorporated  entirely  by  the  application  of  a 

surface- integrated  normal-direction  matrix  addition,  it  was 

shown  that  the  multipolar  moments  for  the  charge  and  current  distributions 
are  related,  which  reflects  the  charge  conservation  on  the  surface. 
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FIGURE  CAPTIONS 

Figure  1.  Geometry  for  the  evaluation  of  the  limit  r+-*r.  Around  r  on 

the  surface,  we  divide  the  surface  into  an  infinitesimal  circle  of  radius  R, 
and  the  rest  of  the  surface.  Then  the  integrals  are  split  up  accordingly. 
The  normal  vector  points  from  the  surface  into  the  dielectric  and  is 
multiplied  by  6>0.  We  denoted  r+-r'  by  a  and  r-r'  by  b.  The  limit  r+-*r 
implies  R>>6>0  and  R-K).  It  appears  that  an  integral  over  the  small  circle 
remains  finite  whenever  the  gradient  of  the  Green's  function  occurs  in  the 
integrand.  This  gives  rise  to  the  first  term  on  the  right-hand  side  of  Eq. 
(3.1). 


Figure  2.  Illustration  of  the  surface  multiplicity.  From  the  origin  0 

in  the  direction  fi,  we  find  points  on  the  surface  which  have 
a  distance  £(0)^,  to  0.  Therefore,  a  description 

of  the  surface  in  spherical  coordinates  requires  a  set  of 
functions 
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